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ABSTRACT 

The  need  for  efficient  and  environmentally  safe  refrigerants  is  of 
great  concern  to  the  Navy,  and  is  driving  them  to  explore  the 
performance  of  alternate  refrigerants  for  use  in  shipboard 
evaporators.  Each  time  a  new  refrigerant  is  considered,  one  must 
also  consider  how  oil  contamination  will  affect  its  performance. 
The  present  study  surveys  the  literature  for  a  predictive  technique 
to  calculate  the  nucleate  boiling  heat  transfer  behavior  of 
refrigerant-oil  mixtures  that  is  fundamental  to  the  operation  of 
refrigerant  systems.  There  were  three  expressions  found  that  could 
be  considered  for  use.  These  expressions  for  the  nucleate  boiling 
heat  transfer  were  empirically  based  formulations,  which  means  that 
a  new  series  of  experiments  must  be  conducted  every  time  a 
different  refrigerant-oil  combination  is  considered  for  use.  An 
expression  that  is  currently  under  development  for  mixtures  in 
general  was  found  that  is  based  only  upon  the  thermophysical 
properties  of  the  fluids  and  their  mixtures.  No  empirical  constants 
are  necessary.  The  mass  diffusion  coefficient  of  the  oil  in  the 
refrigerant  is  found  to  be  the  key  property  to  be  measured  in  order 
to  evaluate  the  predictive  capability  of  this  expression. 
Experiments  are  proposed  to  obtain  this  key  property  and  then  to 
test  the  predictive  capability  of  this  expression.  Modifications  to 
this  expression  will  be  developed  as  the  experiments  indicate  are 
necessary. 
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INTRODUCTION 

The  energy  crisis  of  the  1970 's  brought  an  awareness  of  the  need  to 
design  more  energy  efficient  systems.  The  emergence  of  the 
environmental  concern  over  the  greenhouse  effect  during  the  1980 's 
brought  an  awareness  of  the  need  to  find  and  use  less 
environmentally  toxic  substances  in  thermal  design,  especially  in 
the  design  of  refrigeration  systems.  These  two  forces  have  had  a 
major  impact  on  the  air  conditioning  and  refrigeration  industries, 
and  have  caused  the  Navy  to  search  for  more  energy  efficient,  less 
environmentally  toxic  refrigeration  fluids  for  the  refrigeration 
systems  which  are  on  board  ships  and  submarines. 

Refrigerant  R-114  has  been  the  recent  focus  of  research  at  the 
Naval  Postgraduate  School  [1-4]  because:  (1)  its  operating  pressure 
is  close  to  atmospheric  pressure;  (2)  it  is  highly  stable  with 
respect  to  temperature;  and  (3)  it  is  relatively  low  in  toxicity. 
Under  this  research  program,  extensive  nucleate  boiling  experiments 
have  been  carried  out  to  explore  three  important  topics  relative  to 
the  use  of  R-114  in  large  refrigeration  systems.  First,  the  boiling 
heat  transfer  characteristics  of  pure  R-114  have  been  determined. 
Second,  the  effect  of  oil  contamination  up  to  10%  on  its  heat 
transfer  characteristics  have  been  revealed,  and  third,  the  effects 
of  boiling  surface  enhancement  on  heat  transfer  have  received 
attention.  These  studies  have  provided  considerable  understanding 
of  the  various  phenomena  involved  with  the  increase  in  heat 


transfer  due  to  enhancement  effects  and  with  the  degradation  in 
heat  transfer  due  to  the  introduction  of  contaminant  fluids  in  the 
boiling  process.  There  have  not  been  any  fundamental-based 
predictive  models  developed  that  have  been  used  to  describe  these 
processes  for  various  refrigerant-oil  mixtures,  operating 
conditions  or  enhanced  surfaces. 

Recently,  refrigerant  R-124  has  received  attention  in  the  Navy 
because  it  is  believed  to  be  less  toxic  to  the  environment  than  R- 
114.  Its  operating  pressures  are  reasonable,  and  its  degradation 
properties  relative  to  temperature  are  believed  to  be  acceptable. 
If  R-124  is  to  be  used  in  ship-board  refrigeration  systems,  Zerol 
300  oil,  which  is  different  from  that  used  in  the  previous  studies, 
will  most  likely  be  employed.  Obviously,  if  there  were  a  good, 
fundamental-based  predictive  model  of  nucleate  boiling  heat 
transfer  that  includes  the  effects  of  surface  enhancement  and  oil 
contamination,  it  would  be  of  great  value  as  this  new  refrigerant- 
oil  combination  or  other  alternative  combinations  are  being 
considered. 

Modeling  of  the  process  of  multicomponent  fluids  has  received 
little  attention  as  compared  to  boiling  of  single  pure  fluids.  In 
all  cases,  it  has  been  found  that  the  presence  of  a  second 
component  degrades  the  heat  transfer  as  compared  to  either  of  the 
pure  components,  but  that  modification  of  the  surface  from  the 
normal  "smooth"  character  can  either  increase  or  decrease  the  heat 


transfer  [5].  There  are  no  fundamental-based,  truly  predictive 
expressions  that  have  been  developed  to  describe  the  boiling  heat 
transfer  process  for  multicomponent  refrigerant-oil  fluids,  even 
for  smooth  surfaces  much  less  enhanced  surfaces  [6]. 

The  goals  of  the  present  study  were  to:  (1)  perform  an  exhaustive 
search  of  the  archival  literature  for  studies  of  heat  transfer  due 
to  boiling  of  refrigerant-oil  fluids;  (2)  evaluate  any  predictive 
correlations  discovered;  (3)  idenitify  appropriate  experiments  to 
enable  the  development  of  the  fundamental-based  predictive  model; 
(4)  design  experimental  facilities  to  verify  the  predictive 
capability  of  the  derived  model  for  smooth  surfaces;  and  (5)  design 
experiments  to  study  the  effects  of  surface  enhancement  on  the 
boiling  heat  transfer  of  the  refrigerant-oil  fluid  combinations. 


LITERATURE  SURVEY 

There  have  been  many  publications  dealing  with  the  presentation  of 
data  for  boiling  heat  transfer  of  refrigerant-oil  mixtures.  There 
have  only  been  four  works  that  have  attempted  to  develop 
correlations  to  describe  the  observed  heat  transfer.  In  this 
survey,  these  four  references  will  be  discussed  in  detail,  but  the 
other  studies  will  be  referenced  only  as  needed  to  support  the 
present  work. 


In  1979  Chongrungerong  and  Sauer  [7]  published  the  first 
correlation  of  heat  transfer  for  the  boiling  of  refrigerant/oil 
mixtures.  They  studied  the  case  where  a  prescribed  heat  flux  was 
applied  to  a  single  smooth  tube.  They  included  the  mixture  Prandtl 
number,  the  heat  of  vaporization  of  the  volatile  substance,  the 
effect  of  bubble  density  compared  to  that  of  the  liquid,  the 
applied  heat  flux  and  the  effect  of  pressure  in  the  development  of 
their  correlation.  They  proposed  the  following  expression  for  the 
heat  transfer  for  the  boiling  of  R-ll/oil  mixtures: 


A=0.0523[(      WX>D ;)o.s69(ij^)c.3»pl.6.s(  D  )-o-^  P,y.579 

[\iLbhfg)  kL  0.01588  Pv 

where: 

(J)f  =  Volume  fraction  of  pure  refrigerant 

D  =  Tube  diameter  (m) 
P  =  Pressure  (atm) 

Q/A  =  Applied  heat  flux  (W/  m2    ) 

2 

h   =   Heat  transfer  coefficient  (W/m  ) 
/iL  =  viscosity  of  the  liquid  (g/m-s) 

cL      =   Specific  heat  of  the  liquid  (KJ/g-K) 

kL        =  thermal  conductivity  of  the  liquid  (W/m-K) 


(1) 


Ahf    =  latent  heat  of  vaporization  (W-s/g) 


pLtV  =   density  of  the  liquid  and  vapor  (g/cm  ) 


The  assumptions  behind  this  model  were  that  the  properties  were 
those  of  the  pure  refrigerant  or  the  refrigerant-oil  mixture  as 
appropriate,  the  mixture  was  an  ideal  mixture  (  homogeneous  and  no 
chemical  reactions)  and  the  operation  of  the  boiling  process  was 
steady.  To  develop  the  correlation,  they  performed  experiments  with 
R-ll  as  the  refrigerant,  and  they  used  3GS,  4GS  and  5GS  oils  in 
concentrations  ranging  from  0%  to  10%. 

They  attempted  to  compare  their  equation  with  a  variety  of  data 
that  were  available,  but  they  were  unable  to  do  this  because  there 
were  problems  in  finding  the  appropriate  mixture  properties  that 
were  used  in  the  various  experiments.  They  therefore  selected  some 
typical  properties  and  then  simplified  their  equation  to  the 
expression  given  below: 

£=6.17  (Q/A)°-55<bf2-65P°-2i 

(2) 

Chongrungreong  and  Sauer  [1]  actually  found  that  equation  (2)  did 
a  better  job  of  predicting  the  heat  transfer  for  the  higher  oil 


concentrations  than  the  more  complete  equation  (1)  did.  They  found 
that  equation  (1)  came  within  43%  of  all  of  the  published  data 
where  the  property  information  was  available.  The  worst  agreement 
came  with  R-ll  and  the  highest  oil  concentrations,  approximately 
10%.  Equation  (2)  only  had  a  maximum  error  of  25%  at  the  10%  oil 
concentrations.  The  better  agreement  of  equation  (2)  is  due  in 
part  to  the  selection  of  their  properties.lt  is  important  to  note 
that  they  did  not  model  the  physical  processes  accurately.  For 
example,  the  assumption  of  the  ideal  mixture  does  not  allow  for 
concentration  gradients  to  build  up  around  the  bubbles,  and 
therefore  the  mass  transfer  process  which  moves  oil  away  from  the 
bubble  and  moves  refrigerant  to  the  bubble,  which  can  be 
significant  at  high  oil  concentrations,  is  not  accounted  for  at 
all. 

In  1982  Stephan  and  Mitrovic  [8]  presented  a  nice  discussion  of  the 
physics  involved  in  the  mixture  boiling  process.  They  assumed  that 
the  refrigerant-oil  mixture  could  be  modeled  as  a  binary  mixture. 
They  proceeded  to  explain  what  is  happening  by  referring  to  a  phase 
equilibrium  diagram  as  shown  in  Figure  1.  Plotted  in  this  figure  is 
temperature  versus  the  mole  fraction  of  the  less  volatile  component 
of  the  mixture,  which  in  this  case  is  the  oil.  The  dew-point  line 
and  the  boiling-point  lines  are  drawn.  They  represent  the  locus  of 
states  which  describe  the  mixture  vapor  as  found  in  the  bubble  and 
the  liquid  at  the  initiation  of  the  boiling  process  respectively. 


To  understand  what  physical  phenomena  are  important  to  include,  one 
must  understand  what  happens  during  the  boiling  of  mixtures. 
Consider  a  mixture  that  has  oil  concentration  XA.  The  bubble  starts 
to  form  at  point  A  on  the  diagram,  but  the  temperature  of  this 
mixture  rises  to  point  B  as  a  result  of  the  boiling  process.  This 
creates  a  situation  where  the  liquid  at  the  bubble  interface  has 
concentration  shown  by  point  C  and  the  vapor  is  at  point  D.  In 
reality,  the  formation  and  release  of  the  bubble  from  the  surface 
causes  the  generating  surface  to  cool  slightly,  and  as  a  result  the 
actual  "bubble  state"  moves  to  a  lower  temperature  at  a  slightly 
elevated  oil  concentration  as  seen  at  point  E.  The  surface  then 
heats  up  again  and  the  bubble  returns  to  state  B.  This  process 
repeats  itself  as  the  boiling  process  continues  [8],  When  the  dew- 
point  and  bubble-point  lines  are  close  together,  the  variations  in 
concentration  between  the  vapor  and  the  liquid  at  the  interface  are 
small.  On  the  other  hand,  when  the  curves  are  widely  separated,  the 
effects  of  oil  concentration  gradients  can  become  very  substantial. 
This  is  expected  to  be  the  case  for  most  refrigerant/oil 
combinations  that  are  commonly  used. 

One  assumption  that  is  usually  made  in  modeling  the  mixture  boiling 
process  is  that  the  mixture  concentration  is  uniform  around  the 
bubble,  as  we  saw  in  reference  [7].  Adjacent  to  any  surface  which 
is  transferring  heat  to  its  surroundings  is  a  thermal  boundary 
layer  of  some  finite  thickness.  When  the  bubble  is  of  such  a  size 
as  to  be  significant  relative  to  the  boundary  layer  thickness,  we 
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Figure  1.     Vapor-Liquid  Equilibrium  Diagram 


can  be  sure  that  the  assumption  of  uniform  bubble  temperature,  and 
as  a  result  the  assumption  of  spherical  shape,  are  not  strictly 
correct.  As  a  result,  there  are  nonuniform  flows  around  the  bubble, 
and  the  heat  transfer  around  the  bubble  is  not  uniform.  All  these 
observations  result  in  the  realization  that  the  concentration  of 
the  fluid  surrounding  the  bubble  is  not  uniform  either  and 
therefore  models  that  make  no  attempt  to  consider  the  effects  of 
variations  in  mass  transfer  around  the  bubbles  will  be  doomed  to 
fail  at  higher  oil  concentrations. 

In  1984  Hahne  and  Noworyta  [9]  presented  a  new  correlation  for  heat 
transfer  for  nucleate  boiling  of  refrigerant/oil  mixtures.  They 
developed  an  empirical  expression  that  related  the  heat  transfer 
coefficient  to  the  applied  heat  flux,  and  in  doing  so,  took  the  oil 
concentration  to  be  the  only  variable  other  than  the  applied  heat 
flux.  All  explicit  effects  due  to  temperature  variations,  variable 
properties  and  such  were  contained  in  three  empirical  coefficients 
as  shown  in  Equation  (3). 

h  =  0  .  085  [exp  ibxw)  +exp  (b2w)  ]  qio.B9 -bw) 

(3) 

where  h,  q  and  w,  the  oil  mass  concentration,  are  in  SI  units. 

They  used  refrigerant  R-ll  and  Clavis  G100,  Clavis  G68  and  Oil  22- 
12  in  concentrations  up  to  20%.  In  this  equation  Jb-, ,  b2  and  B  are 
empirical  coefficients  with  units  inverse  to  those  of  the  mass 

10 


concentration.  Hane  and  Noworyta  give  values  of  these  coefficients 
for  each  refrigerant-oil  combination  studied. 

As  one  might  expect,  their  empirical  predictions  agree  quite  well 
with  the  data  since  they  were  not  trying  to  make  the  coefficients 
general  for  all  mixtures.  The  main  drawback  with  this  approach  is 
that  the  experiments  must  be  conducted  with  every  new  refrigerant- 
oil  combination.  It  is  also  only  valid  over  the  concentration  range 
tested;  the  predictions  are  only  as  good  as  the  experiments  where 
the  coefficients  were  determined;  and,  the  only  parameter  beside 
the  applied  boundary  condition  is  the  oil  concentration.  It  leaves 
out  explicit  descriptors  for  the  physical  parameters  that 
characterize  the  physical  phenomena  of  the  boiling  process. 

In  1984  Jensen  and  Jackman  [10]  presented  a  model  of  nucleate 
boiling  of  refrigerant/oil  mixtures  which  does  the  best  job  of 
accounting  for  the  physical  phenomena  that  has  been  published  to 
date.  They  began  with  the  assumption  that  the  variable  property 
effects  were  the  most  important  phenomenon  to  account  for.  They 
assumed  that  there  was  at  the  time  of  bubble  departure  an  oil 
enriched  layer  around  the  entire  bubble.  Underlying  this  assumption 
is  another  assumption,  that  the  refrigerant  is  the  vaporized 
substance  contained  inside  the  bubble.  They  chose  the  time  of 
bubble  departure  as  a  parameter  because  there  is  information  on 
time  to  departure  and  because  at  that  time  the  bubble  is 
approximately  round.  They  assumed  that  the  mass  boundary  layer  was 
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uniform  despite  the  fact  that  the  bubble  was  still  touching  the 
surface  at  the  time  of  departure. 

In  order  to  account  for  the  variation  in  oil  concentration  around 
the  bubble,  they  used  the  concept  of  an  effective  oil  concentration 
that  was  higher  than  the  bulk  concentration.  The  expression  that 
they  used  was 

+  Cpm(6d2z+12dz2+8z3) 
+  pm(6d2z  +  12dz2+8z2) 


(1-C) 
Cpvd: 


(1-C) 


where 


C  =  bulk  mass  oil  concentration  (kgoil/kgmixture) 
pv     =   Density  of  the  vapor  in  the  bubble  (g/  cm2    ) 

d  =  Diameter  of  bubble  at  departure  (m) 

z  =  Thickness  of  mass  boundary  layer  around  bubble  (m) 

This  expression  was  simplified  to 


Ceff=C(l+0.0317AT^a7t53)  (5) 


with  temperature  in  K. 
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This  value  for  the  effective  oil  concentration  was  then  used  in  the 
heat  transfer  expression 


-A  =exp  (  -4  .  095Ceff-55  .  lictff)  (6) 


where  the  value  of  h       in  SI  units  is  given  by  Forster  and  Zuber 


[11] 


*  0.79  0.45  0.049  0.25 

Az  =  0.  00122  1^       Cpl     PJ ^_]A7^afAPs°a7t5  (7) 

a°-^5-29^24p°-24 


As  stated  previously,  this  expression  for  boiling  heat  transfer  of 
refrigerant-oil  mixtures  has  the  strongest  base  in  the  physics  of 
the  boiling  phenomena.  While  this  is  so,  it  is  interesting  to  note 
that,  for  R-113  and  R-ll  refrigerants  with  four  different  oils, 
this  expression  has  an  absolute  error  when  compared  to  the  data 
presented  of  29.6%.  The  authors  attribute  the  error  to  the  scatter 
of  the  data  and  the  lack  of  information  on  mixture  properties, 
especially  the  mass  diffusivity. 

Jensen  and  Jackman  state  that  mass  diffusion  effects  are  very 
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important,  especially  at  higher  oil  concentrations.  They  also 
conclude  that  viscosity  and  surface  tension  are  the  two  properties 
that  they  measured  that  behaved  least  like  the  ideal  mixture  would. 
The  mass  diffusion  process  is  the  least  understood  property  and  is 
key  to  accurate  predictive  models  in  the  future. 

It  is  clear  from  the  above  discussion  that  none  of  these  models 
have  accurately  taken  into  account  all  the  significant  physical 
phenomena  that  are  present  in  nucleate  boiling  of  refrigerant-oil 
mixtures.  The  effects  of  mass  diffusion,  thermal  diffusion,  the 
interactions  of  these  two  processes,  variation  of  mixture 
properties,  transfer  processes  around  the  bubbles  prior  to 
departure  from  the  surface  and  surface  geometry  effects  are  all 
open  questions. 


PHYSICS-BASED  MODEL 

Because  there  have  been  no  models  developed  to  date  for 
refrigerant-oil  mixtures  which  take  the  important  physical 
phenomena  and  parameters  into  account,  it  is  appropriate  to  turn  to 
the  literature  for  general  mixtures.  The  fundamental  physical 
processes  are  the  same  for  all  mixtures.  The  important  thing  is  to 
identify  which  processes  dominate  the  boiling  of  each  combination 
of  fluids.  As  stated  by  Thome  [6],  there  is  no  model  published  to 
date  that  is  comprehensive  in  its  approach  to  modeling  all  the 
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physical  processes  that  go  on  in  the  refrigerant-oil  mixture 
problem. 

The  approach  that  has  included  the  most  physics  in  its  derivation 
for  mixtures  in  general  is  a  study  that  was  developed  by  Bajorek 
[12]  and  is  currently  in  preparation  by  Bajorek,  Lloyd  and  Thome 
[13]  .  This  model  begins  with  the  binary  mixture  models  put  forth  by 
Schlunder  [14]  and  by  Thome  and  Shakir  [15].  First,  it  is  assumed 
that  a  refrigerant/oil  mixture  can  be  considered  to  be  a  binary 
mixture  of  pure  fluids.  This  is  a  reasonable  assumption  since 
refrigerants  and  oils  do  not  react  with  each  other  to  any 
significant  extent  at  temperatures  characteristically  found  in  the 
nucleate  boiling  range.  Assume  that  the  mixture  is  adjacent  to  a 
surface  that  is  hotter  than  the  boiling  temperature  of  the  mixture 
at  that  pressure. 

At  a  nucleation  site  a  bubble  will  begin  to  form.  This  bubble  is 
the  vapor  of  the  refrigerant  since  it  is  significantly  more 
volatile  than  the  oil  under  consideration.  At  the  vapor/liquid 
interface  of  the  bubble,  the  concentration  of  the  refrigerant  is 
depressed  below  that  of  the  bulk  due  to  the  mass  transfer  process 
which  must  take  place  to  move  the  refrigerant  to  the  bubble 
surface,  and  at  the  same  time,  the  concentration  of  the  oil 
increases  above  that  of  the  bulk  as  required  by  continuity.  This 
would  be  uniform  around  the  bubble  if  these  processes  were 
controlled  simply  by  diffusion  of  mass  and  energy.  Because  this 
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takes  place  at  a  surface  and  therefore  in  a  boundary  layer  in  a 
gravitational  field,  convection  becomes  a  player  in  the  boiling 
process,  and  the  temperature  and  species  concentrations  will  vary 
around  the  bubble.  To  properly  account  for  all  these  effects  would 
be  a  very  difficult  analysis.  We  must  simplify  the  analysis. 

Consider  the  schematic  diagram  of  a  bubble  forming  at  a  surface  as 
shown  in  Figure  2.  In  the  actual  problem,  the  bulk  temperature  of 
the  mixture  varies  over  the  diameter  of  the  bubble;  the  heat 
transfer  varies  over  the  surface  of  the  bubble;  the  mass 
concentration  of  the  bulk  fluid  varies  due  to  the  presence  of  the 
wall  and  the  convective  heat  and  mass  transfer  rates  at  the  surface 
of  the  bubble  vary  around  its  circumference  due  to  the  bouyancy 
created  by  the  variations  in  temperature  and  species  concentration. 
One  must  make  some  assumptions  in  order  to  begin  to  solve  this 
problem. 

Following  the  analysis  of  Schlunder  [14]  and  Thome  and  Shakir  [15]  , 
one  defines  an  ideal  heat  transfer  coefficient  by 

hid=q"/{TrTb)  (8) 

and  an  overall  heat  transfer  coefficient  by 

h  =  q"/{Tw-Tb)  (9) 
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Figure  2.    Model  of  Bubble  During  Growth 


Then    it   can   be    shown   that 


hid  Tw-Tb 


It  was  shown  by  Thome  and  Shakir  [6]  that 


=  [1  +  ±1 *]  (10) 


(Ti-Tb)=   {^      (Xb-Xj)  (11) 


where  xi  and  yf  are  from  Figure  1.  They  then  made  the  approximation 

(ri-ri5)=ATJbp[i-exP(-ijf)]  (12) 

where  ATb  is  defined  in  Figure  1,  and  f  is  a  function  of 
thermodynamic  variables  and  accounts  for  the  degradation  in  heat 
transfer  due  to  the  mass  transfer  process.  This  enabled  them  to 
derive  the  following  equation  for  the  nucleate  boiling  heat 
transfer  coefficient 

h     u.hidrdT.  .  ..  <13) 


={l  +  -7r74^^i-*i)  [1-expT]]}"1 


hiH  Q/A     dx 


Bajorek  et  al  [13]  assumed  the  existence  of  an  evaporation 
microlayer  between  the  surface  and  the  bubble.  Above  that  is  a 
relaxation  microlayer  which  reaches  around  the  bubble  to  where  the 
bulk  fluid  exists.  They  further  assumed  that  the  temperature 
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profile  in  the  fluid,  between  the  wall  and  the  bubble,  in  the 
evaporation  microlayer  is  linear.  Given  those  assumptions  they 

derived  an  expression  for   t|i 


kg' 
pAV\Ad(K/Z)12) 


»=.,^   **,..,„  ,o.s  d«) 


where  the  properties  are  in  SI  units. 

Assuming  (yi-xi)  =  (yb-xb)     ,   then  the  expression  for  nucleate  boiling 

heat  transfer  for  the  binary  mixture  is 


={i+^4[^(yb-xb)[i-exp( ~kQlA  )]]}- 


hid  Q/A     dx   Ji>      b  v     pAW^U)1'2 


(15) 


The  predictions  based  on  this  expression  were  compared  to  the  data 
of  Bajorek  [12]  for  acetone/water  and  ethanol/water  mixtures.  This 
expression  was  able  to  predict  the  data  to  9%  on  the  average  with 
a  standard  deviation  of  32%. 

The  question  of  importance  at  this  point  is  whether  or  not  equation 
(15)  is  based  upon  the  correct  physical  phenomena  that  govern  the 
nucleate  boiling  process  for  refrigerant-oil  mixtures.  This  can 
only  be  determined  by  comparing  data  against  the  prediction.  In 
order  to  make  these  comparisons,  the  surface  temperature,  the  bulk 
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mixture  temperature,  the  bulk  mixture  oil  concentration,  the  phase 
equilibrium  curves,  the  pure  component  boiling  temperatures,  the 
bulk  mixture  thermal  conductivity  and  density,  the  bubble  interface 
temperature,  the  thermal  diffusivity  and  the  diffusion  coefficient 
of  oil  in  refrigerant  must  all  be  known.  We  do  not  know  this 
information  for  even  one  of  the  tested  refrigerant-oil  mixtures, 
and  so  some  very  important  experiments  to  find  the  key  mixture 
properties  must  be  designed.  The  most  difficult  property  is  the 
diffusion  coefficient  for  the  oil  in  the  refrigerant  and  the 
thermal  diffusivity  and  conductivity. 

The  important  point  to  realize  at  this  time  is  that  we  have  an 
expression  for  the  nucleate  boiling  heat  transfer  which  contains  no 
empirical  coefficients.  What  remains  is  to  evaluate  its  predictive 
capability  for  refrigerant-oil  mixtures.  To  do  this  we  must  develop 
the  instrumentation  to  measure  the  required  mixture  properties. 


REFRIGERANT-OIL  MIXTURE  PROPERTY  MEASUREMENT 

In  the  derivation  of  equation  (15)  it  was  assumed  that  the  heat 
transfer  to  the  bubble  takes  place  in  the  evaporation  microlayer 
and  therefore  the  relaxation  microlayer  is  comparatively 
unimportant.  This  means  that  the  ideal  heat  transfer  coefficient, 
hid,  is  simply  the  mixture  thermal  conductivity  divided  by  the 
microlayer  thickness.  It  is  possible  to  measure  this  coefficient  if 
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the  bubble  vapor  temperature  can  be  measured.  This  can  in  principle 
be  done  using  interferometry .  Of  course,  then  the  evaporation 
microlayer  thickness  could  also  be  measured.  These  measurements 
would  represent  a  major  contribution  to  the  field.  In  the  absence 
of  these  experiments  one  could  also  experimentally  determine  the 
vapor  equilibrium  curve  for  each  mixture  and  from  that  estimate  the 
bubble  interface  temperature. 

No  matter  what  is  done,  the  diffusion  coefficient  for  the  oil 
through  the  refrigerant  must  be  determined.  The  rest  of  the 
properties  of  interest  needed  to  test  the  heat  transfer 
relationship  of  equation  (15)  for  some  refrigerants  are  available 
in  the  literature. 


PROPOSED  EXPERIMENTS- 

Diffusion  Coefficients 

As  indicated  above,  the  key  property  that  must  be  determined  to 
predict  the  heat  transfer  for  nucleate  boiling  of  refrigerant/oil 
mixtures  is  the  diffusion  coefficient.  The  measurement  of  this 
property  has  not  been  accomplished  for  any  refrigerant/oil  mixture, 
and  so  any  data  will  represent  a  significant  step  forward  in  this 
important  problem. 
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In  order  to  make  these  measurements,  the  technique  outlined  by 
Bidlack  [16]  will  be  used.  In  this  experimental  technique,  a  Mach- 
Zehnder  interferometer  is  used  to  measure  the  variation  of 
refractive  index  of  the  refrigerant  due  to  the  diffusion  of  the  oil 
into  it.  Conceptually  the  experiment  is  quite  simple.  Figure  3 
shows  a  schematic  of  a  test  cell  where  one  creates  a  sharply 
defined  layer  of  oil  beneath  a  layer  of  refrigerant.  There  is  an 
infinite  gradient  at  the  interface  between  the  two  liquids  at  time 
zero,  and  the  liquids  can  be  assumed  infinite  in  extent  relative  to 
the  interface.  As  time  proceeds,  the  oil  diffuses  into  the 
refrigerant,  thereby  changing  its  refractive  index.  If  one  takes  an 
interferogram  at  time  zero  and  then  at  several  time  intervals,  it 
is  possible  to  measure  the  process  of  diffusion  as  a  function  of 
time.  Figure  4  shows  a  schematic  of  a  typical  interferrogram  at 
some  time  t.  With  proper  analysis  of  this  one  can  obtain  the  binary 
diffusion  coefficient  as  follows. 

Consider  the  case  where  the  interface  between  the  refrigerant  and 
the  oil  is  infinitely  sharp  and  is  located  at  x  =  0.  For  x  >  0  the 
fluid  is  refrigerant.  For  x  <  0  the  fluid  is  oil.  At  time  zero, 
the  oil  begins  to  diffuse  into  the  refrigerant,  and  the  refrigerant 
begins  to  diffuse  into  the  oil.  Because  the  diffusion  of  A  into  B 
is  the  same  as  B  into  A,  the  diffusion  front  moves  out  equally  in 
both  directions  from  x  =  0.  Fick's  Second  Law  describes  this 
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Figure  3.   Diagram  of  Diffusion  Cell  (From  Reference  [16]) 
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Figure  4.    Schematic  of  Typical  Interferogram  at 
Time  t  (From  Reference  [16]) 


process  and  is  given  as  follows: 

d2c_    1  dc 
dx2     D12  dt 

with  the  boundary  conditions: 
Case  1: 

For  x>0        x-°°         tzO  c=cx 


t  =  0         c  =  c1  °°>x>0 


x=0  c=  (cx  +  c2) /2  t*0 

Case  2 : 
For  x<0        x--«>         t^O         c=c2 


t  =  0        c=c2        0>x>-oo 


x=0  c={c1  +  c2)  /2  tkO 


(16) 


The  assumptions  for  this  are  that  the  concentration  dependence  of 
D12  is  negligible  and  that  the  diffusion  gradient  has  the  properties 
of  the  normal  distribution  [16]. 

If  one  notes  that  the  concentration  is  directly  proportional  to  the 
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refractive  index,  n,  then  the  solution  of  this  problem  for  any  time 
t  is  as  follows: 


(c-c0)       i      .  r__ —   (n-nQ) 

=  ierf(x/^D^t)=  °  (17) 


(c2-cx)      2       V   iz     (^2-^) 

If   one   defines   the   magnification   of  the    interferogram  to   be  M,    as 
shown    in   Figure   4,    then   one   can   rewrite   Equation    (17)    as 


1  xj  +xk  2 

12    =Trf      erf-^Jll2+erf-U2kzJ)  (18) 

J  J 

In  this  equation  J  is  the  total  number  of  fringes  from  the  top  of 
the  interferrogram  to  the  bottom,  k  and  j  are  the  local  fringe 
numbers  and  the  x1  quantities  are  the  fringe  spacings  as  measured 
on  the  actual  interferrogram.  They  are  related  to  the  true 
dimensions  through  the  magnification  factor  by  x'=  Mx. 

This  is  simplified  by  simply  defining  the  squared  term  to  be 
GRADSQ.  If  one  plots  GRADSQ  vs  time  then  the  SLOPE  of  the  curve 
provides  the  diffusion  coefficient. 

n    SLOPE 

Di2  = T~  (19) 

AM2  V   ; 

The  diffusion  coefficient  defined  this  way  is  assumed  then  to  be 

26 


that  of  the  mixture  at  the  average  concentration  between  the  upper 
and  lower  fluids  in  the  test  section  [16].  The  accuracy  depends 
upon  the  ability  of  the  experimenter  to  create  a  step  change  in 
concentration  at  the  interface  of  the  two  fluids  as  well  as  to 
measure  the  fringe  locations.  Bidlack  claims  that  the  error  is 
within  1%  [16] . 

Thermal  Diffusivity 

The  thermal  diffusivity  of  the  mixture  must  also  be  determined. 
This  can  be  done  using  standard  Differential  Scanning  Calorimetery 
techniques.  In  this  technique,  a  small  sample  of  the  mixture  is 
placed  in  a  special  cup  that  is  is  turn  placed  in  the  heating 
section  of  the  instrument.  At  time  zero,  a  prescribed  heat  flux  is 
applied  and  the  time  response  of  the  fluid  sample  temperature  is 
monitored.  Since  the  sample  cup  is  calibrated,  the  fluid  thermal 
diffusivity  is  the  only  unknown  in  the  problem.  This  procedure  can 
be  repeated  with  samples  of  various  mixture  concentrations  to 
reveal  the  thermal  diffusivity  as  a  function  of  temperature  and 
concentration . 

It  should  be  mentioned  that  once  this  measurement  is  completed,  the 
thermal  conductivity  can  be  obtained  also,  as  can  the  density  and 
specific  heat. 

Vapor-Liquid  Equilibrium  Curves 

It  should  be  noted  that  the  generation  of  the  vapor-equilibrium 
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curves  for  these  mixtures  should  be  easily  accomplished  also.  This 
type  of  measurement  is  commonly  done  in  teaching  laboratories,  and 
therefore  should  not  present  a  great  problem  to  obtain  accurate 

data . 

EXPERIMENTAL  VERIFICATION 

It  is  clear  that  the  final  step  in  the  development  of  this 
predictive  capability  must  receive  some  attention  in  this  report. 
The  predictive  equation  (15)  was  based  on  the  assumption  that  the 
roughness  and  macro  surface  geometry  were  both  sufficient  to 
provide  nucleation  sites.  The  fluid  adjacent  to  the  surface  is 
infinite  in  extent  so  that  the  bubbles  are  not  impeded  in  their 
growth  or  movement  by  neighboring  surfaces.  To  conduct  the 
verifying  experiments  one  must  simulate  this  condition. 

It  is  recommended  that  the  initial  set  of  experiments  be  conducted 
with  single  smooth  tubes  whose  diameter  is  large  compared  to  the 
typical  bubble  diameter.  A  variety  of  oil  concentrations  and 
surface  heat  fluxes  must  be  investigated  in  the  nucleate  boiling 
regime.  When  this  is  complete  to  the  satisfaction  of  the  principal 
investigators,  the  effects  of  surface  geometry  for  the  single  tube 
can  be  investigated,  as  can  the  effect  of  proximity  of  other  tubes 
as  found  in  tube  banks.  A  systemmatic,  step  by  step  set  of 
experiments  can  help  develop  an  understanding  of  the  physical 
phenomena  that  impact  the  heat  transfer  performance  of  the  various 
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refrigerants. 

Finally,  it  should  be  recommended  that  the  choice  of  refrigerant 
for  the  first  experiments  should  be  one  which  has  already  been 
studied  extensively.  This  will  enable  comparison  to  published  data 
faster  than  if  a  new  refrigerant  were  selected  for  the  first 
experiments.  It  is  clear  that  R-114  should  be  among  the  early 
refrigerants  studied  because  of  its  current  interest  to  the  Navy. 
The  recommended  sequence  of  testing  in  order  to  maximize 
understanding  would  be  single,  smooth  tubes  followed  by  enhanced 
surfaces  and  tube  banks. 


CONCLUDING  REMARKS 

The  search  for  new  refrigerants  that  are  environmentally  safer  than 
those  currently  used  and  which  still  perform  well  from  a  thermal 
perspective  is  very  important.  The  problem  of  predicting  nucleate 
boiling  heat  transfer  preformance  of  refrigerant/oil  mixtures  has 
not  been  here-to-fore  addressed  except  through  the  use  of 
experimental  techniques.  This  study  has  addressed  the  problem  of 
finding  a  theoretical  approach  to  this  problem. 

Three  main  studies  were  found  that  developed  empirical  expressions 
for  the  nucleate  boiling  heat  transfer  of  refrigerant/oil  mixtures. 
While  these  studies  present  useful  expressions,  they  are  based  upon 
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experiments  that  must  be  done  for  each  mixture  so  that  one  can  be 
sure  that  the  empirical  coefficients  are  appropriate  to  the 
mixture.  This  need  to  run  extensive  experiments  to  develop  the 
empirical  expression  for  any  new  refrigerant-oil  mixture  clearly 
creates  a  hurdle  in  finding  the  new  refrigerants. 

To  overcome  this  problem,  the  model  being  developed  by  Bajorek  et 
al  [13]  was  discussed.  The  advantage  of  this  model  is  that  it 
depends  only  on  knowing  the  properties  of  the  fluids  of  interest, 
and  there  are  no  empirical  coefficients  that  need  to  be  determined. 
The  main  weaknesses  of  this  model  are  that  it  has  only  been  tried 
out  on  a  very  small  number  of  mixtures,  and  it  does  not  include  the 
effects  of  surface  geometry.  No  refrigerant-oil  mixtures  have  been 
examined  with  the  model,  so  whether  it  properly  accounts  for  the 
physics  of  the  nucleate  boiling  process  of  these  complicated 
mixutres,  even  on  smooth  surfaces,  is  still  an  open  question. 

The  problem  with  trying  this  model  is  that  we  do  not  know  the 
properties  of  the  refrigerant-oil  mixtures.  The  key  property  to 
determine  is  the  mass  diffusion  coefficient,  D12.  An  experiment  to 
determine  this  property  is  proposed.  With  this  capability,  coupled 
with  the  ability  to  measure  the  thermophysical  properties  of  the 
mixtures,  we  should  be  able  to  predict  the  nucleate  boiling  heat 
transfer  for  any  refrigerant-oil  mixture.  Refinements  on  this  model 
also  can  be  proposed  once  the  initial  evaluation  tests  are 
completed. 
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With  an  expression  that  relies  only  on  knowing  the  properties  of 
the  mixture,  we  can  predict  performance  with  sufficient  confidence 
to  make  sound  judgements  on  the  appropriateness  of  new  refrigerants 
for  use  in  Navy  systems. 
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